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Abstract 

A search for the muon catalyzed fusion reaction dd — > 4 He + 7 in the ddfi 
muonic molecule was performed using the experimental [iCF installation TRI- 
TON and NaliTl) detectors for 7-quanta. The high pressure target filled with 
deuterium at temperatures from 85 K to 800 K was exposed to the negative 
muon beam of the JINR phasotron to detect 7-quanta with energy 23.8 MeV. 
The first experimental estimation for the yield of the radiative deuteron capture 
from the ddfi state J=l was obtained at the level n 7 < 2 x 10 -5 per one fusion. 



1 Introduction 

It is understood that investigations of the fusion reactions between hydrogen isotope 
nuclei at low energies are of great importance for determining properties of lightest 
nuclei and for the astrophysics. In particular, there is a need for new or improved 
measurements of many radiation capture reactions included in various astrophysical 
scenarios. Due to the Coulomb repulsion fusion cross-sections cr(E) drop rapidly at 
low (E < lOOkeV) collision energies (in an exponential scale for "bare" nuclei). 
Reaction 

d + d — > 4 Fe + 7 + 23.8Me\/ (1) 
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is involved in both primordial and stellar nucleosynthesis. Its cross section is rather 
small (about 1 pb at 50 keV, compared to 1 mb for main fusion channels d(d,n) 3 He 
and d(d,p) 3 H), and its experimental investigations in dd collisions are rather difficult. 

At energies E > 400 keV reaction (1) proceeds mainly by a d-wave E2 transition 
to the 1 So-state of 4 He J1J. The reason is the identical boson character of the entrance 
channel requiring L + S to be even (L and S are orbital angular momentum and 
total spin of the dd system). At lower energies, the centrifugal barrier suppresses the 
d-wave E2 capture, allowing s-wave E2 transition to the D-state admixture of 4 He. 
Measurements extended to energies below 100 keV have confirmed this picture. 
However, an existence of multipoles other than E2 in reaction (1) was not excluded 
experimentally despite belief that dipole transitions El and Ml with AS = should 
be suppressed due to standard isospin selection rule AT = 0. 

Measurements of the cross section angular distributions a (9), of vector A y and 
tensor A yy analyzing powers performed with a polarized deuteron beam with energy 
E d (\ab) = 80 keV stopping in the target have yielded an unexpected observation of 
the p-wave strength in 2 H(d,^f) A He reaction 0]. It has been found that over 50% of 
the cross section strength at these low energies is due to El and M2 p-wave capture. 
This finding might affect the low energy behavior of the S-function and be considered 
as an isospin selection rule violation. (Some other evidences for non-E2 radiation are 
cited in ||). It would be extremely interesting to observe this p-wave manifestation in 
an independent measurement. 

During past decades experiments in which various fusion reactions between hydro- 
gen isotopes are catalyzed by muons have provided supplementary information about 
these reactions at energies well below the lowest energies accessible by conventional 
beam-target experiments M. In the muon catalyzed (MC) process fusion takes place 
from the bound states of muonic molecules. Nuclei inside muonic molecules are prac- 
tically at rest, being separated by average distances ~ fr 2 /e 2 m^ = 2.5 • 10 -11 cm ( 
m M is the muon mass). 

Muonic molecules can be formed in the states with total angular momenta J = 
and J = 1, that correspond to the relative orbital angular momenta of nuclei L = 
and L = 1. Depending on the hydrogen isotope mixture parameters various states 
of muonic molecules can be populated. This makes possible study fusion reactions at 
super- low energies from prepared s— and p— nuclear states with definite spins. 

Study of the MC fusion process in dd[i muonic molecule resulted in the compli- 
mentary and detailed information about charge-symmetric reactions d(d,n) 3 He and 
d(d,p) 3 H. A significant difference in the p-wave parts of the d(d,p) 3 H and d(d,n) 3 He 
reaction yields was observed in the experiments with low energy polarized deuteron 
beams Comparison of two reaction branches showed some s-wave enhancement 
together with essential p-wave suppression of the proton branch. (This result was then 
interpreted by some as an evidence for charge symmetry breaking forces.) 

Direct measurement of the yields ratio R p (n/p) for reactions 

dd/i — > 3 He + n + fi, ^ 3 He + n, (2) 
dd\i — ► t + p + fi (3) 
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proceeding from the J = 1 state of dd/j, molecules 0, gave the value Rp(n/p) = 
1.42 ±0.03. It agreed with the ratio from |J determined from the elaborate (and model 
dependent) analysis of the in-flight data. Rates of ddfj, fusion reactions (2), (3) from the 
p-wave were experimentally measured || and corresponding nuclear reaction constants 
were extracted from MC data 0. 

The deuteron radiative capture reaction in the cM/i-molecule 

ddfj — > A Hefj + 7 + 23.8 MeV, (4) 

has not previously been investigated because of the extreme smallness of its expected 
yield. In the systematic study of the MC process in deuterium we have recently per- 



formed |L0] measurements in the temperature range T=85-800 K. As in our earlier 



experiments jllj], neutrons from reaction (2) were detected. At temperatures T > 150 



K ddfj molecules are mainly formed in J=l state and fusion reactions proceed from the 
p-wave of relative nuclear motion. Hence, being detected, 23.8 MeV 7-quanta would 
unambiguously indicate a finite p-wave contribution to process (4) rate. 

In view of this, we performed a feasibility study of process (4) detection in our 
last measurements of the <i<i/i-molecule formation rate For this aim one of two 



usually used neutron detectors |TT|]was removed and a gamma detector was installed 
instead. The level of the radiation background in our installation has been measured. 
We present the first experimental estimation for the yield of the radiative deuteron 
capture in the p-wave from the gM/z- molecule. 



2 Experimental setup 

The experimental setup (its layout is shown in Fig. 1) has been described in detail 
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|T2|| . The muon beam of the JINR phasotron was directed into a high pressure 
deuterium target (HPDT) |y|. Scintillation counters 1-4 detected incoming muons. 
Cylinder-shaped counter 5 served to identify the muon stop in the target and to detect 
electrons from muon decay. A coincidence between signals of counters 5 and le, 2e 
served as the electron signal. 



A large neutron detector ND (volume of NE-213 v = 12.5/) [14]], [|T^| was aimed to 



detect neutrons from reaction (2). To reduce a background, the n — 7 separation was 
realized by comparing the signals for the total charge and the fast component charge 
of the ND pulse. The efficiency of the 7-quanta discrimination was better than 10 -3 
for energies E 1>e > 100 keV. 
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Figure 1. Experimental layout. 



The 7-quanta were detected with a NaI(Tl) crystal of a 150 mm diameter and a 100 
mm height. It was calibrated with 7-sources m Co (total energy of two 7-s 2.5 MeV), 
Pu-Be (E 1 = 4.43 MeV) and with 5.5 MeV 7-s from the reaction 

pdfi ^ 3 Hefi + 7. (5) 

Reaction (5) was observed in the test exposures when the target was filled with H/D 
mixture containing about 20% of protium. The energy resolution of the detector was 
15% FWHM at 5.5 MeV. 

Linearity of the energy scale was checked under different voltages supplied to the 
detector in the measurements with available 7-sources 60 Co and Pu-Be. In the used 
amplitude region it proved to be linear at the level 2-3%. Expected position for 7-s 
from (4) is then approximately 200-th channel. Stability of the spectrometric system 
was controlled with 7-sources during the run. 

We estimated the 7-quanta detection efficiency using values of 7-s cross sections in 
Nal and the known solid angle of the detector. Taking into account efficiency losses 
(30 — 40 %) due to the bremsstrahlung in the target walls we determined the efficiency 
for 24 MeV 7-quanta detection 

e 7 = (5 ± 1)%. (6) 

Primary selection of the events detected by the neutron and 7-detectors was realized 
by the trigger. It allowed only those events for further time and amplitude analysis 
which were connected with electron registration, i.e., delayed coincidences \i — n, 7 — e 
were used. Under this condition the timing sequence of the Nal and ND signals was 
registered by flashes ADC and recorded on the PC. 

3 Measurements 

During the run eight exposures were performed at different deuterium temperatures 
and densities. Experimental conditions for them are presented in Table 1. Deuterium 
density is given in relative units: (f> = n/n where n = 4.25 • 10 22 nucl/cm 3 is the liquid 
hydrogen density (LHD). For all exposures the intensity of muons detected by counters 
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1-4 was 2.5 • 10 3 s 1 . The electron counting rate was 15-30 per second depending on 
the deuterium density. 



Table 1. Experimental conditions. 



Exposure 


T, K 


Content, % 


0, 
LHD 


N e 


N n 




C„(H) 


C d (D) 


1 


85 (5) 


20.7 (0.5) 


79.3 (0.5) 


0.84 (0.03) 


712 300 


4 000 


1.2 • 10 5 


2 


110 (5) 


20.7 (0.5) 


79.3 (0.5) 


0.84 (0.03) 


474 600 


4 700 


1.2 • 10 5 


3 


230 (5) 


20.7 (0.5) 


79.3 (0.5) 


0.83 (0.03) 


433 200 


15 000 


4.5 ■ 10 5 


4 


301 (4) 


20.7 (0.5) 


79.3 (0.5) 


0.83 (0.03) 


443 700 


20 200 


6.1 ■ 10 5 


5 


299 (4) 


20.7 (0.5) 


79.3 (0.5) 


0.47 (0.02) 


388 900 


13 900 


4.2 ■ 10 5 


6 


298 (4) 


0.1 (0.1) 


99.9 (0.1) 


0.50 (0.02) 


232 500 


18 100 


5.7- 10 5 


7 


548 (10) 


0.1 (0.1) 


99.9 (0.1) 


0.50 (0.02) 


240 000 


19 500 


5.1 ■ 10 5 


8 


791 (15) 


0.1 (0.1) 


99.9 (0.1) 


0.49 (0.02) 


315 000 


20 500 


6.1 ■ 10 5 



The number of detected electrons from muon decay, N e , was determined from the 
fit of the electron time spectra taking into account the background from muon stops 
in the target walls. The latter one was found in the exposure with the empty target. 
Number of electrons detected per 10 hours of phasotron operation (one exposure) was 
~ (0.5 — 1.0) • 10 6 . Details of the analysis can be found in \T1 . 



The number of neutrons from reaction (2), N n , was obtained from the analysis of 
the time spectra of the events detected by ND and belonging to the neutron region in 
the n — 7 plot |TJ]]. The neutron background was measured in the special exposure with 
the target filled with helium. Thus determined numbers of N e and N n are presented 
in Table 1. 

Exposures made with H/D mixture allowed detection of reaction (5), which was 
used both for energy calibrating and checking the 7-quanta detection efficiency. 

As expected, exposures 1 and 2 are characterized by a low neutron/electron ratio. 
In other words, only small fraction of the muon stops in the target lead to formation 
of dd/x-molecules and subsequent reaction (2), detected in our experiment. It is due 
to the fact |]], [0, ||, ]TI[ that, at low temperatures T < 150 K, the dd\i formation 
rate A^ • < 0.1 • 10 6 s _1 is small compared to the d/x-atom disappearance rate 
^dfi = A + \ddv • ■ (1 — C p ) + Xpd^ ■ 4> ■ C p , where A = 4.55 • 10 5 s" 1 is the free muon 
disappearance rate and = (5.53±0.16) • 10 6 [13] is the pd/j, formation rate. This 
allowed use these exposures for the estimate of the accidental background. Exposures 
3-8 were accepted for the search of 7-s from reaction (4) and estimation of its relative 
yield. 

From the measured numbers of neutrons N n , known efficiency of neutron detection 



e n = 13 % |16j and partial probability of reaction (2) (3 = 0.58 ddp formation 



rates were determined for each exposure. With thus obtained values A^ and data 
from Table 1 we calculated average number of catalysis cycles, n c , using a simplified 
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formula describing /iCF kinetics in H/D mixture (see, e.g. 
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) 



n c = Xddfi • • (1 - C P )/[\ Q + X ddtl ■ 4> ■ (1 - C p ) ■ w dd + Xpdfj, ■ (p-Cp- zupd}. (7) 

Here vdm — 0.07 and w p d = 0.85 are branching ratios of fusion reactions with muon 
sticking to helium with respect to all fusion channels in muonic molecules dd\x and pdfi. 
These reactions lead to the muon loss from the catalysis cycle. Using thus determined 
n c and measured numbers of electrons N e we could calculate numbers of dd\i molecules 
N dd p, formed in each exposure, as = N e ■ n c . The results are presented in Table 1 
and the total number of ddfi molecules for exposures 3-8 was used for the estimation 
of reaction (4) yield. 

4 Analysis of 7-events 

Of all events detected by the 7-detector those with 7 energy 

£ 7 > 17 MeV (8) 

were selected for the further analysis. These events accumulated in exposures 3-8 were 
displayed in a two-dimensional plot 7 time (t 7 ) - electron time (t e ), shown in Fig. 2. 




Figure 2. Two-dimensional t e — 1 7 
plot for the events with E 1 > 17 
MeV detected by Nal detector in 
exposures 3-8. 



2 4^ 6 8 10 12 14 16 18 

t + 

X e , MS 

Fusion events from reaction (4) should arrive after the muon entrance (t^) and 
before the muon decay (t e ), so for the primary selection the following time sequence 
t e , i 7 > tfj,, t e > £ 7 was required, corresponding to the dashed area in Fig. 2. 

It is seen that a noticeable fraction of events in this plot is concentrated at small 
t e — tp, t 1 — tp. These events might be a manifestation of the muon stops in the 
target walls. In their material (Ni, Fe) muon disappears after the average time = 
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0.2 /is, either due to its decay starting the false trigger, or due to the nuclear capture 
with emission of capture products. To reduce the background originating from such 
processes events corresponding to fast 7- and electron emission should be excluded 
from the consideration by introducing a time delay with respect to t^. 

From the other side, time distribution of events resulting from the dd/i molecule 
fusion (4) should obey the exponential law 

/ 7 (t) = Const ■ exp([-X + (j)\ d d^dd]t), 

so allowing a large time delay would lead to the loss of the efficiency 

To provide the suppression of the accidental background and simultaneously avoid 
the efficiency losses the following time intervals were chosen: 

1 jis < t 7 — tfj, < 3 fis; 0.5/is < t 7 — t e < 4.5 {is. (9) 

The corresponding region is indicated by a BOX in Fig. 2 with 3 N* events inside. 

To estimate the background, we selected the area t 7 > 1 fis, 0.5 jis < t e < 4.5 {is 
and found 7 events there. This corresponds to the number of background events in the 
region (9) JV* 1 = 2 ± 1. 

In addition, for an independent estimation of the accidental background, events 
from exposures 1,2 satisfying to (8) and (9) were selected. The number of such events 
normalized to the number of electrons in exposures 3-8 was found N^ 2 = 2. It proved 
to be at the level of the previous estimate obtained from exposures 3-8. 

The background is found to exceed the measured intrinsic background of the in- 
stallation, corresponding to the cosmic ray intensity at sea level, by a factor of 2. We 
conclude that additional background is correlated with phasotron operation. 

The energy distribution of the events detected by the Nal detector and selected 
with criteria (9) for exposures 3-8 (solid line) is shown in Fig. 3. The dashed line is 
the spectrum for the normalized background. 
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Figure 3. Amplitude 7-quanta spectra for the events selected with criteria (8) 
for exposures 3-8 (solid line) and the normalized "background" ones (dashed line). 
Response function of the Nal detector is represented by the gaussian. 
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It is seen from the figure that these spectra practically coincide for energies E 1 > 17 
MeV. Some excess of events for lower energies can be ascribed to the background 
induced by neutrons from reaction (2). 

From considerations above the number of candidate events can be obtained 

iV 7 = N* - iV 7 b = 1 ± 2 (10) 
The measured yield of reaction (4) per cM/z-molecule is evaluated as 

where N*gL is the total number of Nddn molecules accumulated in exposures 3-8: 

Nfl = 3.4 ■ 10 6 (12) 

Using the estimation (6) for the efficiency of 7-quanta registration and taking into 
account the selection efficiency due to the accepted criteria (9) one obtains the detection 
efficiency of 7-s from reaction (4) 

e 7 = (3 ± 0.5) % (13) 

Substituting values (10), (13) and (14) into Eq. (11) we obtain for the absolute 7 
yield per one cM/i-molecule 

rff = (1 ± 2) • 10- 5 . (14) 

In the second run new measurements with deuterium target and NaI(Tl) detector 
of larger size have been conducted and result of the similar analysis of the new data 
set is in agreement with (14) 

?4 2) = (0.8 ± 1.5) ■ 10" 5 . (15) 

Combining (14) and (15) one obtains (at 90% C.L.) 

r? 7 < 2 • 10" 5 . (16) 

From here an upper limit for the radiative fusion rate A 7 from the J=l state of 
dd/u, molecule can be deduced using the experimental value of the total fusion rate 
A} = 4 ■ 10 8 s- 1 § 

A}<8-10 3 s -1 . 



5 Conclusion 

The first attempt has been undertaken to estimate the yield of radiative capture re- 
action (4) from the J=l state of ddfi muonic molecule. The background conditions 
were evaluated and appropriate methods of data analysis were elaborated. The sen- 
sitivity of the present experiment is not enough to make a decisive conclusion about 
the p-wave contribution to the process of radiative dd capture. (With the data from 
H we would expect t/ 7 ~ 10~ 6 .) We plan to achieve this level of sensitivity using 
new 7-detectors of larger efficiency. Of crucial importance is the understanding of the 
background structure and elaboration of background suppression methods. 
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